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Abstract. Geometric characteristics of random fields are exploited to test the 
consistency of density perturbation model spectra with COBE data. These CMB 
maps are analyzed using the number of anisotropy hot spots and their boundary 
curvature. CMB maps which account for instrumental effects and sky coverage are 
Monte Carlo generated. These simulations show that a scale invariant Harrison- 
Zeldovich primordial Gaussian density fluctuation spectrum is consistent with the 
data. The CMB fluctuation coherence angle, based on boundary curvature, gives a 
spectral index n = 1.2 ± 0.3. 

1 Introduction 

Gravitational collapse models of structure formation, predict that at large angular 
scales the cosmic microwave background (CMB) will exhibit random fluctuations. 
The gravitational potential on the surface which last scatter CMB photons (Sachs 
& Wolfe 1967) is responsible for these fluctuations. Thus, the primordial matter 
density spectrum determines the amplitudes of these fluctuations which in k— space is 
P{k) — |<5fc| 2 and can be represented by a power law, P(k) cx k n , with spectral index 
n. The inflationary model requires a scale invariant spectrum with n = 1 (Harrison 
1970; Zeldovich 1972). 

P(k) can be inferred from the CMB temperature autocorrelation function (acf). 
Alternatively, P(k) models can be tested statistically using the hot spots in CMB 
maps. The most important advantage of this second approach is that, unlike acf 
computations, topological analysis does not assume an ergodic radiation field (Cayon 
et al. 1991; Scaramella and Vittorio 1990), and so provides normalization independent 
power spectrum information. This later analysis can also be used to diagonse non- 
random CMB map components (e.g. foreground signals); as an unbias coherence angle 
estimator of an underlying random field (Adler 1981); and to verify the Gaussian 
nature of primordial fluctuations and thus test structure formation models relying 
on primordial seeds (e.g. cosmic strings (Bennett et al. 1992) or texture (Park et 
al. 1991)). Also the existence of non-Gaussian sky temperature contributions from 
non-ergodic behavior can be determined using Monte Carlo simulations. 

The first year COBE-DMR maps (Smoot et al. and references therein) are topolog- 
ically analyzed to check their consistency with the Harrison-Zcldovich power spectrum 
which a preliminary hot spot analysis found to be the case (Gurzadyan and Torres 
1993). 
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2 The geometry of random fields 



The random field excursion set is the domain of all points in which the field takes on 
values T > T„ — va, where a and v are the temperature field standard deviation and 
threshold, respectively. Topological descriptors which have been used to characterize 
CMB excursion regions (Sazhin 1985; Vittorio and Juszkiewics 1987; Bond and Ef- 
stathiou 1987; Coles and Borrow 1987; Coles 1988a; Coles 1988b; Gott et al. 1990; 
Martinez- Gonzalez 1989) are: The number of spots with temperature above T„, and 
below — T v , N v ; the mean area of these spots, A v ; the total spot boundary curvature, 
or genus G„; the contour level lengths, L v \ the total spot excursion area, a„; the 
number of up crossings (one dimension) and the Euler-Poincare characteristic (two 
dimensions (Adler 1981)). 

Excursion set geometry theory allows some descriptors to be derived analytically, 
(e.g. For two-dimensional, stationary, isotropic, random, Gaussian fields the geomet- 
ric properties depend only on the acf and its second derivative at zero lag.) A more 
compact representation is the coherence angle, 2 c — — C(0)/C"(0). The mean number 
of spots on the 47r surface, the mean genus, and the total excursion set area are: 

{ V "<erfc(zVV2)' 
2\ 1/2 v 



{Gv) = {n) ^ ex P(-y)' and ( 2 ) 

(a v ) = 27rerfc(i//\/2). (3) 

The total excursion set area is independent of 9 C and so can be used to identify 
systematic effects in the data. 



3 Computing topological descriptors 

COBE data is apportioned to pixels by projecting the celestial sphere onto a circum- 
scribed sky cube, each cube side being divided into 1024 pixels (Torres et al. 1989). 
The resulting 6144 pixels, each have w 2.9° x 2.9° angular dimensions comparable to 
COBE-DMR angular resolution. The pixel face and Cartesian coordinate identifier 
allow efficient algorithms to be constructed so that the data can be accessed in a form 
close to that of the original. 

Since hot spots are defined in terms of a threshold, u, a binary map is constructed 
for each threshold level, a 1 appearing in pixels whose temperature is greater than or 
equal to T v . Hot spots are located on these bit maps by forming tree data structures. 
When traversing the sky cube the first 1 pixel encountered becomes the root node 
and all contiguous 1 pixels become nodes of that tree. A tree is completed when all 
neighboring pixels are 0. The hot spot area is the number of nodes in the tree. The 
contour level length is found by counting the neighboring 1 pixels of each node. The 
total curvature index, or genus, of a bit map is the number of hot spots minus the 
number of holes (Gott et al. 1992) which are those spots found using the hot spot 
algorithm recursively on the negative of the bit map minus one. 

Simulated Gaussian random field maps arc generated and their hot spot number 
densities and genus, areas, and contour lengths evaluated. Figure 1 shows the simu- 
lated total area a mean values with theoretical a v values while Fig. 2 does the same 
for genus. The agreement is good. 
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Excluding the equatorial band to remove contributions from our galaxy introduces 
a boundary thus enhancing G v and N v at low thresholds. To account for this artifact 
all Monte Carlo generated data sets are treated identically to that of the experimental 
data. 

4 Monte Carlo simulations 

A COBE-DMR CMB map can be simulated by assigning to each pixel a tempera- 
ture given by an expansion of real spherical harmonics as defined in Smoot et al. 
(1991) weighted by the beam filter coefficients We- The finite beam width, cr^, of the 
instrument enters in We as a high-pass filter (Vittorio, Scaramela 1988): 



For random Gaussian fields, the harmonic coefficients of the expansion are random 
variables with zero mean, are Gaussianlly distributed, and have variance determined 
by the primordial matter density perturbation spectrum. 

On the scale of the DMR beam (> 2°) only fluctuations of the gravitational poten- 
tial on the surface of last scattering affect the CMB isotropy and so the Sachs- Wolfe 
effect alone suffices to determine the variances. For power law spectra, P (k) cx k n , 
the Bond-Efstathiou formula (1987) for the variances is adequate since only contribu- 
tions for I up to <~ 30 are important after beam width filtering and further Gaussian 
smoothing. Instrumental noise and sky coverage are accounted for by adding to each 
pixel noise equals to a Gaussian random number with dispersion equal to the noise 
level per measurement divided by y/N b s , where the number of observations is in- 
cluded in the original sky maps. To estimate noise for a single measurement the 
observed temperature squared is apportioned among bins of equal width in 1/N b s , 
and the average in each bin is used. The slope of the (T 2 ) vs. 1/N plot gives the 
noise level squared. The levels for each channel at each of the three DMR frequencies 
were within 7% of nominal values. 

Two maps are Monte Carlo generated, each with its corresponding noise level, 
from which sum and difference maps are constructed and 2.9° Gaussian smoothed. 
Each produced sky map is one realization of the ensamble from which G„, N v , A Vl 
and L v are extracted. After 400 realizations the respective means and dispersions are 
computed. 

Figure 3 shows the means and their 1 sigma error bars for the number of hot 
spots as seen by the 53 GHz DMR radiometer observing a universe with primordial 
Harrison- Zeldovich fluctuations. Also shown are the spot numbers obtained from 
maps made with pure instrumental noise. The separation between these two curves 
indicates both the power of a statistical test and the discrimination efficiency of a 
topological descriptor. Noisy maps have noise and signal curves in close proximity. 
L v and A v have such close noise and signal curves and so can not discriminate be- 
tween power spectra with indices near one. Genus and hot spot number are the only 
descriptors used to test hypothesis. 

To compare COBE-DMR data with cosmological models \ 2 an d the cumulative 
probability distribution of x 2 are constructed from Monte Carlo data (Gott et al. 
1992). For each Monte Carlo map the topological descriptor values U v evaluated at 
thresholds v : 1...25 are used to compute, 




(4) 
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where (Uj) and <jj are calculated from the Monte Carlo simulations for that cosmo- 
logical model. Table 1 shows the x 2 and probabilities for the 53 and 90 sum maps 
and their weighted average. 

5 Preparation of data 

In this analysis four of the six DMR sky maps - 53-A, -B, 90-A, and -B - are used. 
The data is in the sky cube format using Geocentric Ecliptic coordinates. Before 
the topological analysis the dipole and 2fiK kinematic quadrupole components are 
removed, the maps are rotated to galactic equatorial coordinates, and to reduce noise 
the maps have an additional Gaussian smoothing of 2.9°. An equatorial band of 
|b| < 30° is excluded from the analysis. Cosmological hot and cold spots will exists in 
all four maps as unresolved sources, so a weighted average of the four maps is made 
to reduce the noise. 

As a data integrity and analysis software check the diploc, rms 10° sky variation, 
and quadrupole are computed and compared with COBE results. The dipole and 
quadrupole arc obtained by a least squares fit of the harmonic expansion to the map 
and are in agreement with Smoot et al. (1992). 

6 Analysis and Conclusions 

Figure 1, a plot of a„ versus v, illustrates what would be obtained for a pure random 
field independent of its coherence angle. In the same figure, a versus threshold for 
the 31 GHz DMR sum map, shows large deviations which increase with decreasing 
galactic cut indicating galactic contamination at high latitudes. For the 53 and 90 
sum maps no deviations beyond 1 sigma of Fig. 1 are seen above 25 degrees. 

Genus and number of spots were used to test for the consistency of a primordial 
Harrison- Zeldovich perturbation spectrum. The null hypothesis that the structure in 
COBE maps is due solely to instrumental noise is rejected with high confidence as 
seen in Table 1 and Fig. 3. The contribution from known infrared, radio, and x-ray 
point sources is negligible (Bennett et al. 1993) so the observed structure must be of 
cosmological origin. 

Genus and number of spots computed on Monte Carlo realizations of CMB maps 
with a Harrison- Zeldovich power spectrum normalized to 16^iK Q rms is found to be 
consistent with DMR data to a high confidence, as seen in Table 1, as are other 
spectral indices. To see the genus dependence on the spectral index, several sets of 
Monte Carlo realizations are generated with fixed Q rm s normalization and variable n. 
For each index the genus curve as a function of threshold v is fit using (3) to obtain 
the best 9 C . The coherence angle dependence on spectral index is almost linear, as 
seen in Fig. 4. With 9 C errors obtained from Monte Carlo simulations the best value 
of n is found by fitting the COBE genus curves to determine their coherence angle 
which for the 53 GHz sum DMR map is 4.5° ± 0.1° which implies a 1.2 ± 0.3 spectral 
index in agreement with Smoot et al. (1992) using a acf fit. 
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TABLE 1 

X 2 and its probability for number of spots and genus of COBE-DMR maps 



Map 




P(x z > X^) 




p{x z > xU) 


53 


0.62 


95.4 


0.72 


83.3 


90 


1.48 


13.9 


1.73 


8.4 


5390 


1.03 


47.1 


1.03 


47.1 


53* 


6.75 


< 1/400 


7.51 


< 1/400 


90* 


3.59 


< 1/400 


4.38 


< 1/400 



NOTE. - The first three rows are the x 2 statistic (5) for the number of spots xjsj 
and genus Xq an d their probabilities (%) when the Harrison- Zeldovich hypothesis 
is assumed. The 5390 entry refers to the statistics of the weighted averaged map of 
DMR 53 and 90 sum maps. Last two rows give the same statistics when the null 
hypothesis is assumed in the Monte Carlo procedure. 
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FIGURE CAPTIONS 



Figure 1. Mean value of the total area of excursion divided by 27rerfc(V/v / 2) at 

different threshold levels. The horizontal line correspond to the expected 

theoretical value for a random field, independent of its coherence angle. 

Open squares with 1 sigma error bars are for Monte Carlo generated maps. 

Circles are for the 31(A+B) DMR maps with a 30° (filled in circles) and a 

10° (open circles) galactic cut. 
Figure 2. Mean values of the genus as a function of threshold for Monte Carlo 

simulations of a random gaussian maps. The solid curve is the analytical 

function. 

Figure 3. Number of spots normalized to the area of the sphere including cold 
and hot spots. Upper curve is for noise maps alone (open squares). The 
other curve (open squares) is for maps generated with noise and cosmologi- 
cal fluctuations with a Harrison- Zeldovich spectrum (16 [iK normalization). 
Superposed is the data obtained from the 53(A+B) DMR maps (filled in 
triangles). 

Figure 4. Coherence angle versus spectral index as obtained from the Monte Carlo 
runs. Solid curve is a fit to a parabola. 
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